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Abstract : 

In this brief report, we revisit analytical calculation [Mishra, et al, Phys- 
ica A 323 (2003) 453 and Mishra, NewYork Sci. J. 3(1) (2010) 32.] of the 
persistent length of a semiflexible homopolymer chain in the extremely stiff 
chain limit, k — > (where, k is stiffness of the chain) for directed walk lattice 
model in two and three dimensions. Our study for two dimensional (square 
and rectangular) and three dimensional (cubic) lattice case clearly indi- 
cates that the persistent length diverges according to expression (1 — 5c)" 1 , 
where g c is the critical value of step fugacity required for polymerization of 
an infinitely long linear semiflexible homopolymer chain and nature of the 
divergence is independent of the space dimension. This is obviously true 
because in the case of extremely stiff chain limit the polymer chain is a one 
dimensional object and its shape is like a rigid rod. 
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1 Introduction 

The persistent length of a polymer chain measures correlations in the ori- 
entation of the segments of the chain along its length. In other words, the 
persistent length is a measure of a distance along the chain length at which 
the configuration of the chain on an average has memory of the orientation 
of its specific segment. The bending rigidity and thus the persistent length 
is a consequence of short range atomic and molecular interactions present in 
the polymer chain. Since, the persistent length is stemming from the bend- 
ing rigidity of the polymer chain and it can exhibit enormous variation in 
the magnitude. Therefore, if persistent length associated with the polymer 
chain is much smaller than the overall length of the chain, such a chain is 



said to be flexible and stiffness of such chain is unity. When stiffness of the 
chain is approaching to zero, the persistent length of such chain being com- 
parable to it's length and the chain is said to be rigid. However, if stiffness 
of the chain has value in between 0-1, the chain is said to be semi-flexible. 
Actin filaments, microtubules, DNA, protein and collagen are the examples 
of the semiflexible polymers. The persistent length plays an important role 
in describing elastic properties of a semiflexible polymer chain and also plays 
vital role in developing theory of polyelectrolytes solutions. 

Due to excluded volume effect, self avoiding polymer chain has memory 
of it's specific segment and initial bias persists along the walk of the chain 
upto a finite distance (for flexible chains) from initial step of the chain. 
Grassberger [1] initially discussed this problem and showed that the persis- 
tent length of a two dimensional self avoiding flexible polymer chain diverges 
with power law. Later, Redner and Privman [2] suggested that this diver- 
gence is logarithmic. However, through MC studies [3] it has been shown 
that the persistent length could be fitted by power law and by a logarithmic 
function. Eisenberg and Baram [4] demonstrated and confirmed that the 
persistent length of a flexible polymer chain converges to a finite value. The 
situation is different in the case when polymer chain is semiflexible and in 
the extremely stiff chain limit the persistent length of a semiflexible polymer 
chain diverges. 

The aim of present report is to take into account correlations prevailing 
between two distant segments of an extremely rigid polymer chain of an 
infinitely long length in the bulk and to demonstrate through simple calcu- 
lations that the persistent length of such polymer chain when expressed in 
terms of critical value of step fugacity in the extremely stiff chain limit (i. e. 
k — > 0) diverges as a simple pole and the nature of divergence is independent 
of space dimensionality. 

This report is organized as follows: In Sec. 2, we define directed walk 
model in brief and revisit the results of calculation of the persistent length 
for two dimensional (square and rectangular) and three dimensional (cubic) 
lattice to investigate the divergence of the persistent length of an infinitely 
long linear semiflexible homopolymer chain in the extremely stiff chain limit. 
Finally, in Sec. 3, we conclude the discussion by summarizing the results 
obtained. 

2 Model and method of calculations 

We consider following two cases of directedness [5] of the polymer chain for 
square, rectangular and cubic lattices: In the case (i) partially directed self 
avoiding walk (PDSAW) model, the walker is allowed to walk along ±y 
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and +x directions on a square or a rectangular lattice while in the cubic 
lattice case walker is allowed to walk along ±y, +x and +z directions. In 
case (ii) fully directed self avoiding walk (FDSAW) model, the walker is 
allowed to take steps along +x, +y directions in the square and rectangular 
lattice case while along +x, +y and +z directions for the case of a cubic 
lattice. A partially directed self avoiding walk is shown graphically on a two 
dimensional rectangular and a square lattice in figure (1). 




Figure 1: In this figure a partially directed self avoiding walk of a linear 
semiflexible polymer chain is shown (1A) on a square lattice of 9 steps and 
(IB) on a two dimensional rectangular lattice of 11 steps. The step fugacity 
of each step is shown by g, k[= exp(—(3eb)] is the stiffness of the polymer 
chain, ej, is the value of bending energy required to produce one bend in the 
chain and (3[= jr^f] is inverse of thermal energy. The Boltzmann weight of 
the walk shown in figure (1A) is g 9 k 6 and of (IB) is g n k s . 

The partition function of the chain is defined as follows: 

N=oo 

Z(9,k)=J2 £ g N k N » (l) 

7V=0 all walks of N steps 

Where, Nf, is the number of bends in a walk of a polymer chain of N steps 
(monomers) and g is fugacity associated with each step (monomer). The 
partition function of the chain is calculated [6, 7] by us using method of 
generating function technique [5]. 

The persistent length is defined by Mishra el al. [6], as an average length 
of the polymer chain between two successive bends, i.e. l p =< L > / < 
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N b >= ( g dLo g [Z( 9 M y( k dLog[Z(g,k)] ^ where length Q f the chain ig L ( = N ^ 

a being the lattice parameter and N is number of monomers in the chain). 
We have taken value of lattice parameter unity for mathematical sake. 

2.1 PDSAW model on a square lattice: 

The partition function of a linear semiflexible homopolymer chain for this 

model is written as Z PD - S (g,k) = i^g+g^-tg^ ' t 6 ]> wnere 9 is step fu- 
gacity and k is stiffness weight associated with each bend of the polymer 
chain. 

The critical value of step fugacity required for polymerization of an in- 
finitely long linear semiflexible homopolymer chain is determined from the 
singularity of the partition function. The critical value of step fugacity for 
partially directed self avoiding walk model of the chain on a square lattice is 
written in terms of k as, g c = 1+ ^ k [6]. This allows us to write k in terms 
of 5cas , k = ^-. 

The persistent length of the polymer chain for PDSAW model on a 
square lattice can be written as [6], 



2\/2 r r 1, 

:[^2+d ( 2 ) 



P 4 + 3v / 2 L k s 



Substituting k = in Eq. (2), we obtain expression of the persistent 
length as, 



l p = (l-g c )-i (3) 

2.2 FDSAW model on a square lattice: 

For fully directed self avoiding walk model on a square lattice the partition 
function of the chain is written as Zp£,-s(g,k) = i_^ k ~ )g , [6] while g c = 

j-^j:, [6]. Therefore, we have expression for k in terms of g c as, k = 1 ~ 9c , 
while persistent length for this case is, l p = 1 + k~ l , [6]. Substituting the 
value of k in terms of g c for this case too, we get, 

l P = (l-g c y 1 (4) 

2.3 PDSAW model on a two dimensional rectangular lattice: 

We have considered a rectangular lattice which has lattice parameter one 
unit along x— axis and two unit along y— axis. This rectangular lattice can 
be derived from a two dimensional hexagonal lattice and the lattice is shown 
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in figure (IB). The partition function of the polymer chain for this case is 
written as [7]: 

7„ n /_ k \ _ 3g+2g 2 +2g 2 fc-33 + 4g 3 A : -4g3fc2 
^PD-R{9,K) - l-g>-2gW 

In the case of a two dimensional rectangular lattice, the critical value 
of step fugacity for polymerization of an infinitely long linear semiflexible 
homopolymer chain is written in terms of k as, g c = -^==p , [7]. In other 

words, k in terms of g c is written as k = 1 Z 9 ^ , while the persistent length 

^gc 

has dependence on k as, l v = 1 + for PDSAW model on a rectangular 
lattice. The persistent length (on substitution of k in terms of g c ) is re- 
written in terms of g c as, 

l P = (l + 9c)-\l-9c)- 1 (5) 

2.4 FDSAW model on a two dimensional rectangular lattice: 

The partition function of the polymer chain for this case is Zp£)-R{g,k) = 

2g+fl2+ f_ fc ^-^ 3fc " g3fc2 » [ ? ] and we have > 9c = 7=1 , [7] from the singularity 
of the partition function. In this case, k in terms of g c is written as, k = 
X ~ g 9 t and l p = 1 + £j for FDSAW model on a rectangular lattice in two 
dimensions. On substitution of k in terms of g c for FDSAW model on a 
two dimensional rectangular lattice, we get, 

l p = (1 + g c )-\l - g c )-l (6) 

2.5 PDSAW model on a cubic lattice: 

Thr partition function of the polymer chain for partially directed self avoid- 
ing walk model is Z PD - C (g, k) = (1+k lf k 2 } 4 g \ 9 jl^ 2)g+1 , [6]. In this case the 
persistent length of the polymer chain is written as [6] , 

l 2ix 1 [fc- 2 + fc- 1 -4] 

P (l- v / T7 + 2fe- 1 )n 2 + (85 + 21^/17)^-2 

where 

Ul = 85 + 19^17 - (102 + 26 V / T7>- 1 + (34 + 8^/T7)/c~ 2 
and u 2 = 204 + 52VT7 - (272 + 64>/l7)fc -1 . 

The critical value of step fugacity for this case is g c = 2{k+i^-H'^) ' 
For this case too, we follow the method discussed above and substitute, 
k = il-a^Hzn to obtain, 

l P = (l-g c r 1 (8) 
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In this case dependence of the persistent length on k (as shown in Eq. 
(7)) is more involved than the cases discussed in sub-sections (2.1-2.4) and 
expression of the persistent length reduces to a simple form, as we have 
discussed in sub-sections (2.1-2.4), when the persistent length is expressed 
in terms of g c i. e. Eq. (8). 

2.6 FDSAW model on a cubic lattice: 

The partition function of the polymer chain for FDSAW model on a cubic 
lattice is written as ZpD-c(g,k) = 1 „( 1 3 ^ 2 fc)g i M- The critical value of 
step fugacity is g c = njlfc) an< ^ *he persistent length is l p = 1 + ^ [6] 
for FDSAW model on a cubic lattice. In this case too (on substitution of 
k = hi the expression of the persistent length) we obtain, 

l P = (l-g c r 1 (9) 

3 Conclusions 

We have used definition of Mishra et al. [6] to investigate nature of the 
divergence of the persistent length of an infinitely long linear semiflexible 
homopolymer chain in the extremely stiff chain limit i. e. k — > 0. In this 
limit the polymer chain is a one dimensional object and average length of 
the polymer chain between its two successive bends diverges as, (1 — <7 C ) _1 - 
In other words, the persistent length diverges as l p ~ (1 — g c )~ l ~ ^ (where, 
q is an integer) for extremely stiff chain limit. 

When persistent length is expressed in terms of fc, the constant of pro- 
portionality will depend on lattice dimension and model. The constant of 
proportionality will have different value for isotropic model to that of di- 
rected walk model. However, when persistent length is expressed in terms 
of g c , we expect that the nature of the divergence of an average distance 
between two successive bends of the polymer chain will remain same for 
directed and undirected self avoiding walk models and constant of propor- 
tionality will have different value for isotropic (undirected) model than the 
directed walk model. The nature of divergence is identical for partially and 
fully directed walk models of the polymer chain for two and three dimen- 
sional lattices. This is due to fact that in the extremely stiff chain limit the 
polymer chain is a one dimensional object and its shape is like a rigid rod. 

The qualitative nature of variation of the persistent length with stiffness 
of the chain has similar variation for directed and isotropic self avoiding walk 
models in two and three dimensions. However, exact value of the persistent 
length of the chain will depend on space dimensions and type of model 
(directed or isotropic) chosen to enumerate walks of the chain [8] . 
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